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Abstract. It has recently been observed that the normalisation of a one-dimensional 
out-of-equilibrium model, the Asymmetric Exclusion Process (ASEP) with random 
sequential dynamics, is exactly equivalent to the partition function of a two-dimensional 
lattice path model of one-transit walks, or equivalently Dyck paths. This explains 
the applicability of the Lee- Yang theory of partition function zeros to the ASEP 
normalisation. 

In this paper we consider the exact solution of the parallel-update ASEP, a special 
case of the Nagel-Schreckenberg model for traffic flow, in which the ASEP phase 
transitions can be intepreted as jamming transitions, and find that Lee- Yang theory 
still applies. We show that the parallel-update ASEP normalisation can be expressed 
as one of several equivalent two-dimensional lattice path problems involving weighted 
Dyck or Motzkin paths. We introduce the notion of thermodynamic equivalence for 
such paths and show that the robustness of the general form of the ASEP phase diagram 
under various update dynamics is a consequence of this thermodynamic equivalence. 



PACS numbers: 05.40.-a, 05.70.Fh, 02.50.Ey 
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1. Introduction 

It is a fact of equilibrium statistical-mechanical life that if one wishes to establish the 
thermodynamics of a model system it is necessary to calculate its partition function. 
When dealing with nonequilibrium stationary states (for example those that carry a 
current) it is not obvious that the object analogous to the equilibrium partition function, 
namely an appropriately defined normalisation of the stationary distribution, encodes 
the thermodynamics in a similar way. 

By now the Lee- Yang zeros of such a normalisation for a number of nonequilibrium 
steady states have been studied in the complex plane of the microscopically irreversible 
transition rates that define a particular model. In particular, various one-dimensional 
driven diffusive systems [H|2|, reaction-diffusion processes jHJEJEj an d an urn model for 
the separation of sand |E] have been investigated and, in each case, the zeros approach 
the real axis at the appropriate transition points (aspects of this work have been reviewed 
in [Jj). Furthermore, the locus (i.e., impact angle) and density of the zeros near the 
transition point correspond to the order of the phase transition (defined in physical 
terms according to whether there is, e.g., phase coexistence or a diverging correlation 
length) in exactly the same way as for the equilibrium systems long ago considered by 
Lee and Yang [HJ E] ■ 

A striking feature of these various analyses is that the role of the equilibrium 
fugacities in the Lee- Yang approach is taken up by the nonequilibrium transition rates 
present in the model. The aim of this work is to explore this connection more deeply, 
the main result being that one can find equilibrium models with the same phase 
behaviour as the original nonequilibrium model, but within which the transition rates 
are transparently fugacities. 

At a rather abstract level, this phenomenon can be understood in terms of a 
graph-theoretic description of the microscopic process ^01 \H\ In this picture, one 
considers a graph in which each vertex corresponds to a microscopic configuration and 
directed edges to the transitions permitted in an elementary timestep. A representation 
of the steady-state normalisation can then be shown to coincide with the generating 
function of spanning in-trees on this graph Explicitly, if the microscopic transition 
rates (or, in a discrete-time process, transition probabilities) are drawn from the set 
{wi, w 2 , • • • , wn}, the normalisation Z can be expressed as 

Z= Yl n( mi ,...,m N )wrwT 2 ---wT , (1) 

mi,...,mjv 

in which n(mx, . . . , TOjv) counts the number of spanning in-trees on the graph of 
configurations with m 8 edges corresponding to a transition that occurs with rate Wi. 
Thus the rates Wi act as fugacities that control the amount to which various classes 
of spanning in-tree contribute to the overall normalisation. One can see there is the 
possibility for different classes of trees to dominate the normalisation as the fugacities 
that weight the trees are varied. One may conjecture that such changes correspond to 
thermodynamic phase transitions in the microscopic model. 
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The representation of the normalisation JTJ is appealing in that it is uniquely defined 
for any stochastic process with a single steady state and polynomial in the variables tOj. 
In principle, this form can be obtained by summing the principal minors derived from 
the matrix of transition rates, wherein each summand corresponds to the steady-state 
weight of a particular configuration fHJ El In practice, however, evaluation of the 
determinants involved is intractable, and in the cases where nonequilibrium steady-state 
distributions have been obtained via specialised methods, one typically finds polynomials 
that are lower in degree than (|TJ) . We believe that in these cases, the two distributions are 
related by a factor common to all the configurational weights which does not develop any 
nonanalyticities in the thermodynamic limit. In this instance, a Lee- Yang analysis gives 
the same results independently of the form of the normalisation studied and therefore 
we do not worry unduly that the normalisation is generally not uniquely defined. 

A greater difficulty with expressions of the form (0) is that the relationship between 
fugacities and physical observables is unclear. Recently, however, Brak et al. jTT] 
elucidated the connection for the asymmetric simple exclusion process (ASEP) with 
open boundaries, a model that was solved exactly and independently in [T3] and [T3] . 
There, the (reduced) normalisation of the ASEP was related to the partition function 
of an equilibrium surface model in which the two transition rates associated with the 
boundaries turn out to be fugacities conjugate to the densities of surface contacts with 
the horizontal axis from above and below. Thus for this model, there is an interpretation 
of the normalisation from which the thermodynamics can be extracted using standard 
equilibrium statistical mechanical techniques. 

This result raises the question of whether the relationship between nonequilibrium 
transition rates and equilibrium fugacities extends to other nonequilibrium steady states. 
In this work we investigate this question by revisiting another exactly- solved variant of 
the ASEP, namely that with parallel lattice updates ^Hl El QE|- Whereas the case 
with a random-sequential updating scheme, referred to above, was introduced to model 
the kinetics of biopolymerisation JH], the parallel version arises as a special case of 
the Nagel-Schreckenberg model of traffic flow (2Hj. From this point of view, the phase 
transitions in the model's steady state can be related to traffic jamming. 

We begin by recalling the definition of the parallel-update ASEP along with some 
key results. In addition to the two boundary parameters of the random-sequential 
variant, there is a further quantity, p, which relates to the degree of parallelism in 
the dynamics. This gives rise to a normalisation with a considerably more complicated 
structure than that for the random- sequential ASEP. It is therefore appropriate to check 
that the Lee- Yang zeros of this normalisation correctly reproduce the known phase 
behaviour for the model in the complex plane of the defining transition probabilities. 
This being the case, we move on to establish partition functions for equilibrium 
surfaces or, equivalently, random walks on two different lattices that have the same 
thermodynamic phase behaviour as the parallel-update ASEP (and thence jamming 
in its guise as a traffic model). In each case the extra parameter, p, generalises the 
surfaces or walks obtained for the random-sequential ASEP in a physically meaningful 
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Figure 1. Dynamics of the ASEP. The arrow labels indicate the probabilities with 
which the corresponding transitions occur; site labels are also indicated. Note that the 
moves are conducted in parallel. 



way. In one case, p is a fugacity related to the density of horizontal segments of walks 
on a triangular lattice and in the other it is related to the density of peaks on the 
surface. In order to obtain these results we find that it is of great benefit to consider a 
grand-canonical ensemble in which the system size is a quantity exhibiting equilibrium 
fluctuations about its mean. By considering the general equilibrium thermodynamic 
properties of the surfaces, we are able to gain a further insight into the robustness of 
the phase diagram for the ASEP under different microscopic updating schemes. 

2. The Parallel-Update ASEP 

The microscopic dynamics of the asymmetric simple exclusion process with parallel 
dynamics are as shown schematically in Fig. 1 and defined as follows. Particles are 
introduced with probability a in each timestep at the start of a one-dimensional lattice 
with L sites (if the first site is empty) and leave with probability /3 at the other end. 
They can hop with probability p one unit to the right if that space is empty (or remain 
still with probability 1 — p) , otherwise if the space to the right is blocked they remain 
stationary. The update is applied to all particles simultaneously in contrast to the 
random sequential update scheme. Nevertheless, the latter is recovered in the limit 
p oc dt — > (where dt is the size of each timestep) if a and f3 are rescaled to remain 
proportional to p. 

Quantities such as the current and density correlation functions for the ASEP 
with fully parallel dynamics (p = 1) were first worked out by Tilstra and Ernst 
Although a cluster approximation was used, the results obtained were believed to be 
correct in the thermodynamic limit (L — ► oo). This claim was subsequently confirmed 
in the independent works of and ^H] which employed two different matrix product 
ansatze to yield exact results for general p. Furthermore, these calculations 

provide the representations of the normalisation required for the present work. 

These arise by finding the steady state weights f(C) for each configuration C that 
satisfy the stationarity condition 

J2 [f{C')W{C> -> C) - f{C)W{C - C')\ = (2) 

in which W(C — > C) is the probability of making the transition from configuration C to 
C in a single timestep. A normalisation is then obtained through 

£ = (3) 
c 
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As remarked in the introduction, the weights are defined only up to an overall scale. In 
the calculations of 17J, a reduced form of the full expression ((TJ was given in terms of 
the quantities 

n \ / n — r — 1 \ / n + 1 
r + t ) [ t j ~~ 



Q>n,r ^ ] 
t=0 



along with the conventions 



1 and 



r + t + 1 



n — r — 2 
t - 1 



With these definitions the normalisation is 
Z L = Z' L + pZ' L _ x , 

where 

(j9(i-/5)//3r +i -(p(i-«)/«r +i 



for any integer X. 



(4) 



I 

r=0 



(5) 



(6) 



(p(l-/3)//3)-(p(l-a0/a) 

This should be contrasted with the simpler expression for the ASEP with random 
sequential updates 

l-r)\(l/(3) r+1 - (l/a) r+1 



j^r(2L 

r=l 



(7) 



L\{L-r)\ (1//3) - (1/c) ' 

in which a and f3 are transition rates rather than probabilities. 

It is also possible to define grand-canonical partition functions ^3 El 121] by 
summing over system sizes with fugacity z in both cases, so that Z p (z) = J2l ZlZ L 
in the parallel case and Z{z) = J2l Zlz l in the sequential one. These definitions yield 



Z r , 



where 



a/3{l+pz) 2(1 -p){a(3-p 2 z) - a(3b 2 (l - pz) - a(3b 2 J(l + pz) 2 - Az 



P 



and 



a(3 

a(p — a) 
p 2 (l — a) 

in the parallel case and 
Z(z) 



2p\l - f3)(l - a)(z - z hd )(z - z ld ) 



[(1 -P) -(!-«)(! - 



Zld 



z hd 



4a{3 



(9) 



(10) 



(12) 



(1 - 2a - VI -4z)(l -2(3- VT^Iz) 

in the sequential case [2U 122 121] ■ 

The phase diagram can be extracted by considering the large L behaviour of Zl and 
Zl, either directly or from the asymptotics of the grand-canonical generating functions. 
In the parallel case using the latter approach [T7j we can see that the leading singularities 
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Figure 2. Phase diagram of the ASEP with parallel-update dynamics. The second- 
order lines are at a, (3 = 1 — ^/l — p and the first-order line in bold a — f3 runs from the 
origin to meet these. Also indicated as a dotted line is the locus 1 — p = (1 — a)(l — (3) 
along which the mean field solution is exact. 



in (jHJ) come from the poles at z id) Zhd m the low- and high-density phases respectively, 
or from the square root singularity, 
2-p-2^T=p~ 



p2 



(13) 



in the maximal current phase. 

The phase diagrams for the random sequential and parallel updates have similar 
topologies. That for the parallel case is shown in Fig. 121 where the main features to 
remark upon are the second-order transition lines at a or (3 = 1 — ^/l — p and the first- 
order line, depicted in bold, a = (3 < 1 — \/l — p. In addition, the mean field solution is 
exact along the dotted locus 1 — p — (1 — ct)(l — (3). The region in which both a and (3 
are greater than 1 — ^/l — p is a maximal current phase, that with a < 1 — a/1 — p and 
(3 > a is a low-density phase, and that with j3 < 1 — ^/l — p and a > (3 is a high-density 
phase. There are additional transitions within these phases in which the bulk remains 
the same but the boundary profiles change and these occur at a or (3 = 1 — ^1 — p. 
In a similar vein the phase diagram for the random sequential updates ASEP shows a 
maximal current phase for a and (3 > 1/2 separated from high- and low-density phases 
by second order transition lines a — 1/2, /3 > 1/2 and (3 = 1/2, a > 1/2. There, the 
high- and low-density phases are separated by a first-order line a = (3 < 1/2. 

3. Lee- Yang Zeros 

It has been observed that, for non-equilibrium systems, a possible equivalent of the 
reduced free energy might be given by 

/ = - lim jZ L . (14) 

L^oo ±j 

Since for parallel-update dynamics the singularities of (jHJ) determine the asymptotic 
behaviour of the coefficients Z^, we can see that these will simply be given by 



/ = In z mc , for a, (3 > 1 - Jl - p 
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f=\nzi d , for (3 > a, a<l-Jl-p 



f=\nz hd , fora>/3, /3<l-yjl-p. (15) 
For random sequential updates the singularities in (JT^jl give 

/= lni, for a, /3 > 1/2 , 

/ = lna(l — a) , for (3 > a, a < 1/2 , 

f = \n (3(1 - (3) , for a > 0, (3 < 1/2 , (16) 

where / is similarly defined as — lim^oo Z L /L. 

It has been suggested that for the random sequential update ASEP and other similar 
non-equilibrium models the Lee- Yang picture of equilibrium phase transitions [HI EJ 
might still apply pQ. The starting point of Lee and Yang's work was the consideration 
of how the non-analyticity characteristic of a phase transition appear given that the 
partition function is a polynomial on finite lattices or graphs. The form of the reduced 
free energy makes it clear that any non-analyticities are associated with zeros of the 
partition function when the appropriate fugacities (e.g. y = exp(— 2h) for a spin model 
in field) are extended into the complex plane. In general, for a lattice with L sites, the 
partition function can be written as a polynomial in the fugacity, 

Z L = J2DrV\ (17) 

r 

where the degree of the polynomial is proportional to L. As the polynomial may be 
completely expressed in terms of its zeros, y r , so too may the (reduced) finite-size free 
energy: 

f L {h) ~ -j)nY[{y - yr{h)) . (18) 

In the thermodynamic limit, L — > oo, the zeros move in to pinch the real axis at the 
transition point (or points) and the entire phase structure of the model is determined 
by their limiting locus and density. From p8j) . one has 

f(h)~-f dyp(y)\n(y-y(C)), (19) 
Jc 

where C represents some set of curves, or possibly regions, in the complex y-plane on 
which the zeros have support and p(y) is the density of the zeros there. 

The expressions for the reduced free energy for the ASEP allow a direct 
determination of the locus of partition function zeros. In general, the partition function 
zeros delineate the boundaries of different phases since they give the loci along which 
the free energy is singular, signalling a phase transition. In the Lee- Yang approach 
one considers an extension to complex parameters, so the free energy, /, becomes a 
complex quantity, fa, which is the metastable free energy per unit volume in the various 
phases given by I = 1, 2.., with Re/; = / characterising the free energy when phase I is 
stable. The loci of zeros are then determined by demanding Ref = Refj for adjacent 
phases Since we have / = In z mCt hd id for the various phases of the parallel-update 
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ASEP, we can immediately see that the locus of zeros for the first-order transition will 
be determined by 

\zhd\ = \zid\ , a,P<l - \jl-p ■ (20) 
Substituting in for Zhd, %id from (fTU|) and (fTTj) . this gives 
g(g ~ P) 

Similarly, the loci of zeros for the second order transitions are determined by \zhd\ = \z, 
and \zia\ = \z mc \ or, from (jTUjl . (|TT|) and (|13)1 

p 2 (l-/3) 
a{p — a) 
p 2 (l — a) 

These formulae are analogous to those observed in the case of the first-order transition 
line in random sequential dynamics 

|a(l-a)| = (23) 

and for the second order lines 

| 7 (l- 7 )| = l/4, (24) 

where 7 = a, (3 depending on the region of the phase diagram. 

For definiteness we now consider the zeros of the normalisation for parallel- 
update dynamics given in (JZj) with a convenient rational value of p = 16/25. This 
places the horizontal and vertical second order transition lines of Fig. 2 at a and 
(3=1 — y/1 — p = 2/5, respectively. If we consider the former, we can choose a fixed 
(3 value, say (3 = 3/5, which will take us across this line as a is varied. Substituting 
p = 16/25, (3 = 3/5 into the second of (J22J gives the analytic locus of zeros in the 
complex a-plane which is plotted in Fig. 3. We have also used (JBJ) to calculate the zeros 
numerically for system sizes up to L = 250, and these are plotted in Fig. 3, showing quite 
good agreement with the limiting analytical results. Similarly we can take a different 
value of (3 which crosses the first-order line, say (3 = 1/10, and use (|2*U|) determine the 
analytic locus, which is shown in Fig. 4 along with the numerically determined zeros for 
a system of size L = 250. 

Referring to the discussion of the partition function zeros for the random sequential 
update ASEP in [2], we thus see the Lee- Yang zeros technique still captures the nature 
of phase transitions with parallel updates. In Fig. 3 the zeros approach the phase 
transition point at a = 2/5 at a angle of tt/4 as might be expected for a second order 
phase transition. In the first-order case in Fig. 4 the zeros approach the transition point 
at a = [3 = 1/10 at an angle tt/2, which is characteristic of a first-order transition. 

It is also possible to look at zeros of the normalisation in the complex p-plane 
for given a and (3. In this case one would expect to see a transition at the larger 
of 1 — (1 — a) 2 and 1 — (1 — (3) 2 . The locus of zeros can be obtained by using, 



a(p — a) 



p 2 (l — a) 



(21) 



2-p-2y/r^ 

p 2 



P 



(22) 
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Figure 3. The locus of zeros in the complex a-plane for the parallel-update ASEP 
with p = 16/25, (3 — 3/5, characteristic of a second-order transition. The diamonds 
show for comparison numerically determined zeros for a system of size L = 250. 



o.i 




-O.i 



Figure 4. The locus of zeros in the complex a-plane for the parallel-update ASEP 
with p — 16/25, (3 — 1/10 characteristic of a first-order transition. Again numerically 
obtained zeros for a system of size L = 250 are superimposed. 

e.g. \(5{p - (3)/(p 2 (l - (3))\ = |(2 — p — 2^1 - p)/p 2 \ for fixed (3 and complex p. The 
resulting cardioid locus is shown in Fig. with the transition appearing as the cusp on 
the real axis. 
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Figure 5. The locus of zeros in the complex p-plane for the parallel-update 
ASEP with (3 = 2/5. The transition point is given by the cusp on the real axis at 
1 - (1 -/3) 2 = 16/25 = 0.64. 



4. Thermodynamically Equivalent Equilibrium Models 

The Lee- Yang analysis of the nonequilibrium normalisation for the parallel- up date 
ASEP in the preceding section suggests that the transition probabilities defining the 
model are related to equilibrium fugacities. To establish this connection concretely, 
we seek equilibrium models which are thermodynamically equivalent to the ASEP with 
parallel dynamics, which we define to occur when 

(i) the extensive part of the equilibrium free energy is identical to that of the logarithm 
of the nonequilibrium normalisation specified by equations flU)-®; and 

(ii) the parameters a, (3 and p are conjugate to physical densities in the equilibrium 



A consequence of this definition is that (given the foregoing Lee- Yang analysis 
of the nonequilibrium normalisation) such equilibrium models will have the same 
thermodynamic phase behaviour. An example of a thermodynamic equivalence was 
recently provided for the random-sequential limit p — > by Brak et al. [HJ EZj . 

Our method of choice in establishing such equivalences in the more general case 
of parallel dynamics is to relate the grand-canonical normalisation (JHJ) to generating 
functions for various types of lattice walks. In an earlier work [21] we demonstrated 
that in the limit p — > this approach is technically simpler than transforming the 
transfer matrices that build up the lattice-path partition functions to the matrix product 
expressions used to solve the ASEP, which was the method used in [TTJ [2JJ . 

To this end we first note that (JHJ) can be expressed as 



model. 




(25) 




1 2 3 4 5 6 



Figure 6. An excursion on the rotated square lattice. Note that between the initial 
and final steps, the excursion comprises sub-excursions away from the line y = 1 and 
that the length of the path along the x-axis is equal to the number of up-down pairs. 



where 



x.{z) -|( I />.:)- v 'Vl. 4 pi) 2 - I: j (26) 



is one of the roots of the pole terms, 
x(p — x) 



(27) 



p 2 (l — x) 

appearing in the grand-canonical partition function. 

Let us begin by briefly recapitulating the limit p — > for which the equilibrium 
equivalence has been established [TT] I2*H I2*7j . since the methods we employ will be of 
utility in the general case. It is clear from (|2l)Jl that the function X- = X-jp remains 
finite in the limit p — > 0. Therefore, to keep (|25jl finite we write 

z (z) = [a/p]\P/p](l-p[x-(z)/p]*) 
A ) ([a/p)-[x4z)/p))Wp]-[x-(z)/p]) ' 1 ) 

which indicates the limit p — > must be taken with a/p = a, (3/p — j3 as claimed in 
Sectional Now note, from (|2Tjl . that the function X-(z) satisfies 

z = x-{l-xJ), (29) 

in the limit p — > 0. This expression implies that x_ is the generating function of the 
number of excursions on the 45° rotated square lattice, i.e. paths that touch the x- 
axis only at the start and end. The fugacity z is conjugate to the length of the path, 
measured in terms of the number of pairs of up- and down-steps — see Fig. El 

The generating function for excursions, Ge(z), can be derived by observing from 
Fig. EH that each comprises an initial up-step and final down-step (which together 
contribute a fugacity z) between which any non-negative number of excursions away 
from the line y = 1 can occur. Therefore Ge{z) can be expanded as 

G E (z) =z(l + G E {z) + [G E (z)} 2 + [G E {z)f + ■■•)= ^— - . (30) 

v ' f — ^eKz) 

Comparison with (|29J) indicates that x = G E ( claimed. 

To complete the path interpretation of (|25J) in the random-sequential limit, we 

recast in terms of the rescaled quantities to obtain 

1 




Figure 7. A one-transit walk obtained by composing two Dyck paths, one above the 
axis and one below. Each contact with the x-axis from above, apart from the start, 
gives a factor 1/a and each contact from below, apart from the end, gives a factor of 
1/13. 



in agreement with |2~T] I24| . Expansion of the denominator reveals that this expression 
describes an ensemble of paths with two sets of excursions, one with a fugacity 1/a per 
excursion (or, equivalently, return to the x-axis) and the other with fugacity 1/(3. For 
convenience we take the first set of excursions to go above the x-axis and the second 
below thus defining a one-transit walk, a realisation of which is shown in Fig. [7| In 
the mathematical literature walks that are constrained to lie on or above the x-axis are 
often called Dyck paths (see e.g. |28j). A one-transit walk is then a composition of two 
Dyck paths. 

We now use similar methods to establish equilibrium statistical mechanical 
interpretations for the parallel- up date ASEP with general p. 



4-1. Walks on the triangular lattice 

Let us now consider excursions on a triangular lattice, often called Motzkin paths (such 
as that shown in Fig. |HJ) in which each horizontal step contributes a fugacity u and 
an up-down pair v. The generating function Gt{u,v) for these excursions can be 
derived as follows (but see also, e.g., [2H])- As before, we note that each excursion 
comprises an extremal up-down pair (contributing a factor v) and in between there may 
be any combination of horizontal steps (each contributing a factor u) and excursions 
(contributing Gt(u,v)), including the null path. That is 

G T (u, v) = v(l+ [u+G T (u, «)] + [u+G T (u, v)f + [u+G T (u, v)} 3 + . . .) 

= V (32) 

l-u-G T (u,v) ' 1 ; 

Solving this equation one finds 



1 - u - J(l - u) 2 - Av 
G T (u, v) = ^ 1 ) (33) 

in which the negative root is taken to ensure all coefficients appearing in the power-series 
expansion in u and v are positive. 

A direct comparison of with (|2fi|) yields an interpretation of the normalisation 
for the parallel-update ASEP in terms of triangular-lattice walks with a fugacity \fz for 



Dyck Paths, Motzkin Paths and Traffic Jams 



13 



y 
1 

1 2 3 4 5 6 * 

Figure 8. An excursion on the triangular square lattice. In this case both horizontal 
segments and sub-excursions occur on the line y = 1. The total length of the path is 
the sum of the numbers of up-down pairs and horizontal segments. 





Figure 9. A one- transit walk on the triangular lattice with contact weights for the 
axis. The horizontal segments have a weight —pz whereas up-down pairs contribute 
a weight z, so that the overall power of z (not shown) is the horizontal length of the 
path. Meanwhile, there are contact fugacities of p/a from above and p/(3 from below. 



each up- or down-step and — pz for each horizontal step. This reads 

Z (z) = l-pG T (-pz,z) 2 f v 

p{ ) (l-(p/a)G T (-pz,z))(l-(p/(3)G T (-pz,z))- { ] 

Note that the power of z in the generating function Gt{—PZ, z) measures the horizontal 
length of the path as shown in Fig. |S] since each up-down pair and horizontal segment 
has the same length in the x-direction. 

We remark that only the denominator is relevant in terms of the leading 
thermodynamic behaviour of the model. One way to see this is to note that the the 
factors in the numerator do not alter the location or nature of the singularities in 
the complex plane. More precisely, if one calculates the canonical (fixed system-size) 
free energy of the walk, it transpires that the numerator supplies only subextensive 
contributions (this point is demonstrated more explicitly in Section 03 below). Hence we 
have established a thermodynamic equivalence, as defined at the start of this section, 
between the parallel-update ASEP and the one-transit walk on the triangular lattice 
with fugacities p/a and pj (5 associated with contacts from above and below, a fugacity 
z conjugate to the path length and a negative fugacity — p conjugate to the density of 
horizontal segments as illustrated in Fig. El 

The origin of the negative fugacity can be understood from the discussion of the 
Introduction. Recall that one form of the normalisation is the generating function 
of spanning in-trees on the graph of configurations with edges weighted according to 
the elementary transition probabilities, i.e. Eq. (pQ). Since this model has discrete- 
time dynamics both the probabilities a,(3,p that particles move and their complements 
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1 — a,l — (3, 1 — p that particles remain stationary appear as weights. We therefore do 
not require an expansion of the normalisation in powers p to have positive coefficients; 
however, we should expect to find a series expansion in p and q = 1 — p with positive 
coefficients. 

One can show this indeed the case for the function Gt{~PZ, z) and thence the 
expansion of the denominator of ()34|) . Expanding (}3*3*j) and making the change of variable 
p = 1 — q one finds after some routine algebra that 

* ( -*"-( i+ gsK-)(-"-iH- (35) 

in which the coefficients of z n q m are clearly positive. The coefficients in (|35|) are the 
Narayana numbers 



N(n,m) = -[ '")[ •"], (36) 




n \m J \m — 1 

which appear in many combinatorial contexts. Of particular note here is the fact that 
they are known to count Dyck paths with 2n steps and m peaks, so this is strongly 
suggestive that the parallel-update ASEP also admits an interpretation involving 
standard Dyck paths with weighted peaks. We shall see that this is indeed the case 
in the next section. 

It is also possible to find a triangular lattice path interpretation of the normalisation 
()25|) in which the quantities p and q appear more naturally as positive fugacities. To do 
this, one observes that, since p + q = 1, (|26|) can be rewritten as 



(1 -pz) - -pz) 2 -4qz 
x^/p — pz= - . (37) 

Now the right-hand side of this expression is equal to GT{pz,qz), i.e. the generating 
function of excursions on the triangular lattice with a fugacity q associated with each 
up-down pair and p with each horizontal step. Expressing (J34|) in terms of Gt{pz, qz) 
yields 

z = 1 + (2p/q)G T (pz,qz) - (p/q)G T (pz,qz) 2 
'" (1 - w a G T {pz,qz))(l - w b G T (pz,qz)) 

with 

«. = , (39) 

qa 

As before, the terms in the numerator of (|38|) are thermodynamically unimportant. 
Therefore the phase behaviour of the parallel-update ASEP corresponds to that of a 
one-transit walk on the triangular lattice with contact fugacities w a and w b from above 
and below respectively, as shown in Fig. E3 

In both of the triangular lattice pictures one can understand the limit p — ► 
intuitively. In this limit, paths with horizontal segments are excluded from the partition 
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Figure 10. An alternative one-transit walk on the triangular lattice. The horizontal 
segments have a weight pz whereas up-down pairs contribute a weight qz. Again the 
overall (omitted) power of z is the horizontal length of the path. The contact fugacities 
from above and below are modified to w a and Wb respectively. 



function and one recovers the one-transit walks on the rotated square lattice as discussed 
above. 

4-2. Weighted-peak walks on the rotated square lattice 

Another natural extension of the one-transit walk on the rotated square lattice that 
arises in the case of parallel-update dynamics has q = 1 — p conjugate to the number 
of peaks along the walk. By "peak" we mean here a maximum above the x-axis, or a 
minimum below. 

We begin, as before, by constructing the generating function Gp(z, q) of excursions 
on the rotated square lattice with the addition of a fugacity, q, associated with each 
peak (a similar derivation is given in |31j). In this instance one has between the initial 
up-step and final down-step (which contributes a factor z) either the null path (creating 
a peak with fugacity q) or any positive number of excursions in which the number of 
peaks is the sum of those from each excursion. Thus we obtain the expansion 

G P (z, q) = z(q + G P (z, q) + [G P (z, q)f + [G P (z, q)} 3 + • • •) (41) 
Solving this equation and discarding the unphysical root, we find 



„. . (1 -pz) - J(l -pz) 2 -Aqz 

Gp(z, q) = , (43) 

which is none other than the right-hand side of (|37j). Hence we have immediately a 
representation of the parallel- up date ASEP normalisation 

zn= l + (2p/q)G P (z,q)-(p/q)G P (z,qy 
P{ ' (l-w a G P (z,q))(l-w b G P (z,q)) ' 1 J 

whose denominator describes a one-transit walk on the rotated-square lattice with each 
peak weighted by a fugacity q and contact fugacities w a and Wb associated with contacts 
from above and below and as given in the previous section as indicated in Fig. ^2 

Again, the numerator of this expression is thermodynamically unimportant and 
we see that the limit p — > removes the weights associated with the peaks. This 




Figure 11. A weighted-peak one-transit walk on the rotated square lattice. Here 
there is an additional weight q associated with each maximum above the cc-axis and 
minimum below. The contact fugacities from above and below are the same as in 
Fig. 033 



corresponds, as it should, to the one-transit walks described at the start of this section 
in the context of random- sequential dynamics. 

4-3. Other update dynamics 

Other update dynamics for the ASEP have been extensively investigated, including 
sublattice parallel dynamics, and backward and forward ordered sequential dynamics 
[SHI IM|- It was observed in [13 HE] that there was a close relation between sublattice 
parallel dynamics, backward and forward ordered sequential dynamics and the (totally) 
parallel dynamics we have discussed so far. For all these the rules for making a hop 
are the same, but in the sublattice parallel case all site pairs i, i + 1 with even i are 
updated on even time steps and all pairs with odd i on odd time steps, whereas in 
ordered sequential updates the sites are updated one at a time starting at the left hand 
side and proceeding to the right (forward ordered) or right to left (backward ordered). 
The canonical normalisation for the sublattice parallel (sp), backward ordered sequential 
(bos) and forward ordered sequential dynamics (fos) is given by the Z L of ©, which in 
turn gives the grand canonical normalisation Z sPt h OS) f os 



a(3 [2(1 - p) (a(3 - p 2 z) - af3b 2 (l - pz) - a/3b 2 ^J(l + pz) 2 - Az } 
Z SP , bos , fos(z) = 2p\l - P)(l ~ *)(z - z hd )(z - z ld ) - ' (45) 

for all three updates [T7] . 

We can again rewrite (J45)) in a form which is recognizably that of a generating 
function for one-transit walks. Using x_(z) — | ((1 + pz) — \J (1 + pz) 2 — 4z^ we have 

aP(l-x-(z)) 
(a — x_(z))(j3 — X-(z)) 

which should be compared to the expression for parallel updates in (J2HJ). With pUj) at 
hand any of the walk representations previously employed may be used to recast the 
grand-canonical normalisation in the form of a one-transit walk generating function. 
For instance, using the triangular lattice generating function one can write 

7 (A (l- P G T (-pz,z)) 

SP ' hOS ' /o4 J (l-(p/a)G T (-p Zj z))(l-(p/P)G T (-pz,z)) ■ 



Aip, bos, fos\Z) — — TZWTa 73\ ' v 4C v 
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Since the singularities in (J47)) arise from the square root in Gt{— pz, z) and the poles 
in the denominator, which are identical to the parallel case, we find the same phase 
diagram. 

5. Thermodynamics of One- Transit Walks 

To finish we briefly investigate the thermodynamics of one-transit walks. It is, in fact, 
possible to derive qualitatively the phase diagram with only a little knowledge of the 
detailed structure of the walks. This provides a picture within which one can understand 
clearly the robustness of the phase diagram of the ASEP under a variety of updating 
schemes [321 EBl 03 ■ 

Consider a walk of length L that comprises n a excursions above and n& excursions 
below the x-axis. For added generality we associate with each walk a number, m, 
that counts some property that depends on the shape of the walk, but not on whether 
a particular excursion lies above or below the axis. The number of peaks is such a 
property. If there is an energy e a (e&) associated with each of the n a (rib) contacts from 
above (below) and an energy e' with the property counted by m, the free energy of the 
walk is 

F = n a e a + n b e b + me' - S (L, n a , n b , m) , (48) 

in which S(L, n a , n b , m) is the entropy (logarithm of the number of realisations) of paths 
with fixed n a , n b and m. 

The consequence of m being independent of whether an excursion lies above or 
below the axis is that 

S (L, n a , n b ,m) = In Q(L, n a + n b , m) , (49) 

where fl(L, n, m) is the number of walks of length L comprising n excursions and fixed 
m. It then follows that if (|48)l can be minimised with respect to n a , n b and m in such a 
way that n = n a + n b is nonzero, this minimum is achieved by all the excursions being 
on the side of the x-axis for which the contact energy is smaller. That is, 

F = n min{e a , e;,} + me' — In Q(L, n, m) . (50) 

If e a = e b , the free energy is neutral to the flipping of an excursion, and so along this 
coexistence line, the crossing point will wander along the x-axis. This is suggestive of a 
first-order phase transition. 

It is intuitively clear that the number of possible walks increases as the number 
of contacts decreases. This implies that once the smaller contact energy has increased 
above some critical value, the entropy gained from avoiding the axis (i.e. becoming 
desorbed) is not offset by any energy saving in contacting the axis. Furthermore, since 
the free energy depends only on the smaller of contact energies, the critical value must 
be independent of the larger contact energy. The effect of e' will be to favour a particular 
shape of path, and thus we would expect this to affect the critical contact energy at 
which the path becomes desorbed. These considerations give rise to a phase diagram 
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Figure 12. Generic phase diagram of a one-transit walk. The density of contacts 
is nonzero on the side of the axis for which the contact energy e aj & is lower. The 
(negative) critical energy e c at which the entropy of the walk causes it to desorb from 
the axis depends on e' and how it is coupled to the shape of the path. 



with the topology shown in Fig. El which is expected be common to a large class of 
one-transit walks. 

For concreteness, let us now consider the specific case of the peak-weighted walk 
illustrated in Fig. [TT] which has m as the number of peaks, e a ^ = — lnu> afe and e' = — In q. 
The number of walks of length L comprising n excursions and m peaks is the coefficient 
of z L q m in [Gp(z, q)] n . This quantity can be obtained readily by applying the Lagrange 
inversion formula to the functional relation (pf2*)) for Gp j^UEHl- This reveals that for 
these walks, Q(L,n,m) is the coefficient of u L ~ n q m in 

n ( u 
7 9 + 



u, 

from which one quickly finds 

{ ^(^\(L-n-l\ <n<m<L 

Q(L,n,m) = < 1 n = m = L (51) 

[ otherwise. 
In the thermodynamic limit, L — > oo, one finds the extensive part of the entropy to be 

s(p,a) = lim — In Q(L, pL, aV) 

L— >oo ±j 

= - [a\na + 2(1 - a) ln(l - a) + (a - p) ln(a - p) - (1 - p) ln(l - p)\ , (52) 

in which p = n/L and a = m/L are the thermodynamic densities of contacts and peaks. 
Minimising the free energy per unit length 

/ = pmin{e a ,e fe } + ae - s(p, a) , (53) 
with respect to these densities one finds their equilibrium values are 
* = q{w-lf-l 
9 (w-l)(q[w-l] + l) ' 1 ] 

a* = ^ W ~ 1} , (55) 
q{w - 1) + 1 ' y J 
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for those values of q and w = max{w a , Wb} where the physical requirement < p* < 
a* < 1 holds. One can check that the inequality p* < a* < 1 is automatically 
true whenever p* > is satisfied. The latter condition, however, is violated when 
w c < 1 + g -1 / 2 . Then the free energy is minimised at some point on the boundary of 
the physical part of the (p, cr)-plane. It is easily verified that this point has p* = (i.e. 
total desorption of the path) and a peak density 

a l/2 

"* = T+g^ ' < 86 > 
One then finds that in all phases the free energy takes the form 

/ = 21n(l-CT*)-ln(l-p*) . (57) 

Using these results together with the definitions (ffilj) and (jlOj) . one recovers the 
expressions for the ASEP free energy given in Section El 

This provides explicit confirmation that the peak-weighted one-transit walk and 
the parallel- update ASEP are thermodynamically equivalent, despite the fact the grand- 
canonical partition function of the latter does not agree exactly with that for the ASEP 
(}44|) which has additional factors in its numerator. The reason for the equivalence 
is that these factors describe a small, fixed number of additional excursions which 
contribute only to the subextensive part of the entropy, and hence become irrelevant in 
the thermodynamic limit. 



6. Conclusions 



The Lee- Yang theory of partition function zeros has been shown to apply to the 
ASEP with parallel- up date dynamics, correctly pinpointing the first- and second-order 
transition lines in the model in the complex plane of transition probabilities. This fact 
suggests that there is an interpretation of these transition rates as equilibrium fugacities. 
Indeed, we were able to make the correspondence concrete by mapping the normalisation 
of the ASEP onto partition functions for equilibrium path problems. We found that for 
parallel updates there were several possible path transcriptions, including two types of 
adsorbing paths on a triangular lattice and a class of peak-weighted paths on the rotated 
square lattice. In all cases the appropriate scaling limit recovered the one-transit walk 
associated with the random sequential update ASEP. 

The benefit of finding such mappings is twofold. First one sees quite clearly 
why a Lee- Yang analysis of a nonequilibrium normalisation should be appropriate 
in the complex plane of transition rates or probabilities: under the mapping to an 
equilibrium model, one finds that they are fugacities. Secondly, standard techniques 
from equilibrium statistical mechanics can be used to learn about the thermodynamics 
of nonequilibrium models. For example, we showed that the structure of the phase 
diagram for one-transit walks is essentially unchanged when an additional fugacity 
associated with the shape of the excursions is introduced. In turn, this provides an 
explanation as to why parallelising the dynamics of the ASEP does not unduly change 
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its phase behaviour. This phenomenon can also be understood in terms of an extremal 
current principle for driven diffusive processes jHUEHl- This is suggestive that that the 
equilibrium walk picture explored in the present work is applicable to a wider range of 
processes to those studied so far. Conversely, other nonequilibrium systems are known 
to be equivalent to equilibrium models (such as the "raise and peel" interface model [37] ) 
so one might expect Lee- Yang theory to apply to these also. Meanwhile, we note that the 
path equivalences have also been utilised to calculate density fluctuations in the ASEP 
through a suitable interpretation of the corresponding matrix-product expressions [38J. 

In this work, we focussed on equilibrium models that are thermodynamically 
equivalent to the ASEP with parallel- up date dynamics. At the start of Section 0] we 
defined this equivalence to exist when the extensive parts of the free energies (defined 
for the nonequilibrium case as the logarithm of the normalisation) coincide for the two 
models. In practice, this meant that the equilibrium surface models we described had 
partition functions that differed from the normalisation of the ASEP. This is reflected 
in, e.g., nontrivial numerator factors in the grand-canonical normalisations (|H4jl . ([«?%]) 
and (jUj) . These factors did not concern us here, since we were interested only in the 
thermodynamic phase behaviour of the models: if one wishes to consider finite-sized 
systems, one should take care to include them. Nevertheless, it turns out that even for 
finite-sized systems, the loci of zeros of the normalisation obtained from an "artificial" 
grand-canonical generating function constructed from the denominator alone are similar 
to those obtained using the full expressions. 

Thus far, our analyses (and related calculations in [TTJ EH]) have relied on the 
prior existence of solutions for the nonequilibrium steady states under consideration. 
This is, of course, an unsatisfactory state of affairs, and ideally we would like to be 
able to establish the path equivalences directly from the definition of the microscopic 
process. We suggest, admittedly speculatively, that the extremely compact forms of 
the grand-canonical normalisations for the ASEP suggest that such an enterprise may 
be possible (though no small challenge). In particular, such a study might yield a new 
way to determine nonequilibrium steady-state distributions. As a first step towards 
this goal, it would be of interest to investigate other nonequilibrium states for which 
thermodynamically equivalent equilibrium models have not yet been derived. Chief 
among these is the partially asymmetric exclusion process (PASEP) jSHJ ED] in which 
particles may hop to the left as well as to the right across the lattice. When the left 
and right hop rates are equal, there is an additional transition to a phase in which there 
is a zero current in the thermodynamic limit. If one could understand this transition 
from an equilibrium statistical mechanical point of view, it would give further weight to 
the growing evidence that phase transitions in equilibrium and nonequilibrium steady 
states are two sides of same coin. 
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